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ABSTRACT 

Aims. Violent gravitational interactions can change the morphologies of galaxies and, by means of merging, transform them into 
elliptical galaxies. We aim to investigate how they affect the evolution of galactic magnetic fields. 

Methods. We selected 16 systems of interacting galaxies with available VLA archive radio data at 4.86 and 1.4GHz and compared 
their radio emission and estimated magnetic field strengths with their star-forming activity, far-infrared emission, and the stage of 
tidal interaction. 

Results. The estimated mean of total magnetic field strength for our sample of interacting galaxies is 14 + 5 yuG, which is larger than 
for the non-interacting objects. The field regularity (of 0.27 ± 0.09) is lower than in typical spirals and indicates enhanced production 
of random magnetic fields in the interacting objects. We find a general evolution of magnetic fields: for weak interactions the strength 
of magnetic field is almost constant (10 - 15 ^G) as interaction advances, then it increases up to 2x, peaks at the nuclear coalescence 
(25yuG), and decreases again, down to 5 - 6fiG, for the post-merger remnants. The main production of magnetic fields in colliding 
galaxies thus terminates somewhere close to the nuclear coalescence, after which magnetic field diffuses. The magnetic field strength 
for whole galaxies is weakly afi'ected by the star formation rate (SFR), while the dependence is higher for galactic centres. We 
show that the morphological distortions visible in the radio total and polarized emission do not depend statistically on the global or 
local SFRs, while they do increase (especially in the polarization) with the advance of interaction. The constructed radio-far-infrared 
relations for interacting and non-interacting galaxies display a similar balance between the generation of cosmic rays, magnetic fields, 
and the production of the thermal energy and dust radiation. 

Conclusions. The regular magnetic fields are much more sensitive to morphological distortions induced by tidal interactions than 
are the random fields. As a result the polarized emission could be yet another indicator of an ongoing merging process. The found 
evolution of magnetic field with advancing interaction would definitely imply a stronger effect of magnetic fields on the galaxy 
surroundings in the earlier cosmological epochs. The process of strong gravitational interactions can efficiently magnetize the merger's 
surroundings, having a similar magnetizing effect on intergalactic medium as supernova explosions or galactic winds. If interacting 
galaxies generate some ultra-high energy cosmic rays (UHECRs), the disk or magnetized outflows can deflect them (up to 23°), and 
make an association of the observed UHECRs with the sites of their origin very uncertain. 
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1. Introduction 

Both observational evidence and theoretical considerations show 
that galaxy interactions are the most important processes dur- 
ing the formation and evolution of galaxies (Struck 1999). This 
is the case both locally in the nearby universe and, in particu- 
lar, cosmologically at high redshifts, where galaxy encounters 
are expected to be more frequent. During such encounters, di- 
rect collisions between stars are highly unlikely, while collisions 
between the interstellar medium (ISM) of the two disks are in- 
evitable and may trigger star formation. The interaction-induced 
starbursts, nuclear activity, feedback, gas infalls, and outflows, 
all play critical roles in shaping ISM and the surrounding inter- 
galactic medium (IGM) (Smith et al. 1201 Oi l. 

Some low-speed collisions lead to merging of galaxies. Gas- 
rich mergers could reach the phase of a luminous infrared galaxy 
(LIRG), or even of an ultraluminous infrared galaxy (ULIRG) 
(Sanders & Mirabel 1996). Violent galaxy collisions are also 
suggested as a source of large-scale shocks and one of pos- 
sible sources of ultra high energy cosmic rays (UHECR; e.g. 
Bhattacharjee & Sigl l2000T) . However, not all galaxy interactions 
will result in a merger, and it is equally important to study differ- 
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ent types of galaxy interaction to understand their role in galaxy 
evolution. 

There were several systematic attempts to study the effects 
of collisions on structures of galaxies and the properties of dif- 
ferent ISM phases. An early work of Toomre & Toomre ( 119721 ) 
showed that numerical simulations quite readily reproduce the 
observed morphology of a wide range of interacting galaxies and 
explain the formation of bars, tidal arms, and bridges. Toomre 
(1977.) compiled a sample of 1 1 ongoing merger systems at var- 
ious stages of interaction, which became known as the Toomre 
sequence. The sequence itself and numerical simulations of its 
members explain how elliptical galaxies can be formed by the 
coalescence of spiral galaxies. The Toomre sequence was stud- 
ied over the years in great detail in various wavelength bands 
in order to understand the evolving properties of interacting sys- 
tems. It was shown that in moving along the sequence from early 
to late-stage mergers, a larger fraction of atomic gas appears out- 
side the optical bodies (Hibbard & Gorkom [l996l ). In the final 
stages of nuclear coalescence, there is little H i within the rem- 
nant bodies, while the tidal material flows in and is able to trig- 
ger a nuclear starburst. In accordance with that. Hummel et al. 
( 1990) found that the central sources of interacting galaxies are 
by a factor of five times brighter in the radio wavelengths than 
those in the isolated spirals. 
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The studies of broadband (V, I) and narrowband (Ha) images 
obtained with the Hubble Space Telescope (HST) by Laine et al. 
(2003 ) show surprisingly little evidence of any distinct trends in 
nuclear properties along the Toomre merger sequence. There is 
just a rise in the nuclear luminosity density in the most-evolved 
members of the sequence. However, the near-infrared (NIR) ob- 
servations with the HST (Rossa et al. 2007), which are less af- 
fected by dust obscuration show a distinct trend for the nuclei 
in the sequence to become more luminous as the merger process 
advances. In a different way than the non-interacting objects, the 
Toomre-sequence galaxies also reveal newly formed stars that 
are more concentrated toward their centres. Rossa et al. also ar- 
gue that if left to evolve and fade for several Gyrs, the merger 
remnants would have the K-hand luminosity profiles of normal 
elliptical galaxies. 

More complicated evolution during the merging process was 
identified by Brassington et al. (2007) in the X-rays observed 
with Chandra. They studied a sample of nine interacting and 
post-merger objects. The original Toomre sequence was ex- 
tended in this work to include very old (up to a 3 Gyr) merger- 
remnants. First, they find that unlike Ha and optical emission, 
the X-ray luminosity peaks about 300 Myr before the nuclear co- 
alescence. The subsequent drop in X-ray emission is likely due 
to large-scale diffuse outflows coming out of the galactic disks 
and reducing the hot gas density. Second, about 1 Gyr after the 
coalescence the remnants are X-ray-fainter than typical elliptical 
galaxies, which are known to be generally strong X-ray emitters. 
Thirdly, the dynamically older mergers seem to rebuild the hot 
galactic halo and rise again the X-ray luminosity. This can be ex- 
plained by an outflowing wind driven into the halo by a declining 
population of Type la supernova (Brassington et al. 2007). 

Although magnetic field constitute an important ingredient 
of the ISM, little is known about how galaxy interactions affect 
its strength and evolution, and how they respond to the evolu- 
tion of other ISM phases presented above . The most detailed 
analysis of magnetic field in disks of interacting galaxies was 
done by Chyzy & Beck (I2004I I for the Antennae galaxies. They 
show that the mean total magnetic field strength (about 2Qfj.G) 
is significantly greater than in isolated grand-design spirals and 
may be due to the interaction-induced star formation. The field 
regularity (the ratio of regular to random field components) ap- 
peared to be lower than in undistorted galaxies, which can indi- 
cate enhanced turbulent motions in the merger system. A strong 
regular component (of lOjuG) was found to trace gas-shearing 
motions along the tidal tail. When and how the enhanced mag- 
netic field arose and when it would disperse in this system are 
still not known. 

The evolution of magnetic field in merging systems like the 
Antennae may look different from the evolution of the X-ray 
component revealed by Brassington et al. (I2007I I. In general, the 
early-type galaxies are typically radio weak, and strong radio 
emission is only observed in ellipticals with an AGN activity. 
In such active galaxies, generation of strong magnetic fields, al- 
though possible in accretion disks, is not connected with any 
large-scale or small-scale galactic dynamo processes, which op- 
erate in typical spiral or irregular galaxies (Beck et al. 119961 
Brandenburg & Subramanian 2005 , Chyzy et al. 1201 11 1. The first 
process requires a strong differential rotation with uniform sign 
of turbulence twisting, the second depends on the ISM turbu- 
lence, which likely declines in early-type galaxies owing to a 
weaker star formation. 

The goal of this paper is to investigate the evolution of mag- 
netic field during galaxy interactions. As magnetic fields are best 
studied from the synchrotron emission, we restrict our consider- 



ations to objects with available radio data. The polarized prop- 
erties of radio emission also provides a unique diagnostic tool 
for recognizing gas flows and their distortions. In Sect. |2] we 
describe the selection criteria of the constructed sample of in- 
teracting galaxies and present the used (VLAQ) radio data and 
the details of their reduction. The next section provides a brief 
description of each interacting system and its properties in radio 
total and polarized emission. We present estimations of equipar- 
tition magnetic field strength for all the galaxies and compare 
them with those for isolated spirals (Sect. [J!9l l. 

The statistical analysis of all the galaxies begins from inves- 
tigating the evolution of magnetic fields as the interaction and 
merging process advance ( 14.11 ). In Sect. 14.21 we examine the re- 
lation of magnetic fields to the star formation rate (SFR). An 
analysis of observed morphological asymmetries in total and po- 
larized radio emission is provided in Sect. 14.31 The radio-far- 
infrared (FIR) diagram for interacting galaxies is presented in 
Sect. 14.41 and compared to the literature results. We also discuss 
the possible magnetization of IGM by magnetic fields of merger 
origin and the propagation effects of UHECRs in the vicinity 
of colliding galaxies (Sect 14. 5t . We summarize our findings in 
Sect.H 

2. The sample and radio data reduction 

To explore the evolution of magnetic field in interacting galax- 
ies, we have compiled a sample of suitable objects. First, we 
chose objects of the best-known merging galaxy sequence of 
Toomre (1 19771 ) for which appropriate radio data are available. 
We restricted our search to the archive of the VLA, frequencies 
4.86 and 1.4 GHz, and D, C, or BC configurations of the array 
to maximize the sensitivity for extended radio emission. There 
are 5 out of 1 1 merging systems with relevant radio data (see 
Table[Tli. Next, we selected six out of nine objects with available 
radio data from the X-ray merging sequence of Brassington et 
al. (2007) (Table [U. Furthermore, we performed a more general 
search for angularly-large interacting objects with the usable ra- 
dio data. We looked for galaxies to fulfil the following criteria: 



- are RC3 galaxies (de Vaucouleurs et al. 1991 ) recognized as 
interacting objects according to Surace et al. (2004 ); 

- are detected in the infrared (e.g. are included in the IRAS 
Bright Galaxies Sample: Sanders et al. 120031 ): 

- at least one object in the interacting system has an angular 
size larger than 2' to study distribution of the radio emission, 
and are also smaller than 9' to be within the primary beam of 
the VLA antenna; 

- declination > -40° ; 

- total flux of at least one galaxy in the interacting system is 
larger than 80mJy at 1.4 GHz (from the NVSS, Condon et 
al. [T998l) . 

We obtained 24 objects satisfying those requirements. For 
ten objects, which are members of seven interacting galaxy pairs, 
we found high-quality VLA archive radio polarimetric data at 
4.86 GHz (see Table[T]). We enriched this sample by a few impor- 
tant types of interacting objects: the polar-ring galaxy NGC 660, 
the weakly interacting spiral NGC 4254 in the Virgo Cluster (see 
Chyzy et al. 2007bl Chyzy |2008), and UGC813/UGC816 (see 
Sect. 13.61 ). In total, in our final sample we assembled 24 pre- 
and post-merger galaxies from 16 recognized interacting sys- 
tems. Although they by no means constitute a complete sample 



' The VLA of the NRAO is operated by Associated Universities, Inc., 
under cooperative agreement with the NSF 
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Table 1. Interacting systems and their members. 



Name 



System name 
or other name 



Type 



Dist. 

[Mpc] 



Inkl. 
[deg.] 



Pos. angle 
[deg] 



H I extent 
[kpc] 



Interaction 
Stage 



NGC 876" 


NGC 876/877 


SAc: sp 


50.8 


77.9 


27.3 


N/A 


NGC 877*' 


NGC 876/877 


SAB(rs)bc LIRG 


50.8 


35.6 


138.0 




NGC 4254 


The Virgo Cluster 


SA(s)c,LINER,HII 


17 


32.0 


60.0 


40 


NGC 2207"' 


NGC 2207/IC 2163 


SAB(rs)bc pec 


35.0 


58.2 


115.9 


56 


IC2163" 


NGC 2207/IC 2163 


SB(rs)c pec 


35.0 


78.5 


102.6 




NGC 5426 


Arp 271, NGC 5426/5427 


SA(s)c pec 


26.7 


69.7 


0.5 


43 


NGC 5427^ 


NGC 5426/5427 


SA(s)c pec,Sy2,HII 


26.7 


25.5 


135.0 




NGC 6907^ 


NGC 6907/6908 


SB(s)bc 


44.5 


37.5 


66.7 


78 


NGC 6908 


NGC 6907/6908 


S 


44.5 


N/A 


N/A 




UGC 12914-5 


Taffy galaxy 


SAB(rs)c,Sbrst,Sy2 


59.6 


54.1 


159.6 


52 


UGC 12915 


Taffy galaxy 


Sdm 


59.6 


72.9 


135.4 




UGC 813 


Taffy2 galaxy 


Sb 


67 


72.0 


110.3 


73 


UGC 816 


Taffy2 galaxy 


Sc 


67 


62.0 


170.0 




NGC 660 




SB(s)a pec,HII 


12.3 


78.8 


11.9 


47 


NGC 403 8* '^^ 


Arp 244,The Antennae 


SB(s)m pec 


26.8 


51.9 


133.2 


133 


NGC 4039* '^^ 


The antennae 


SA(s)m,pec,LINER 


26.8 


71.2 


132.0 




NGC 662 F 


Arp 81 


Sb pec,LIRG 


86.4 


70.8 


142.5 


63 


NGC6622'' 


Arp81 


Sa 


86.4 


27.3 


117.0 




NGC 520"'^ 


Arp 157 


pec 


30.2 


77.4 


130.0 


103 


NGC 3256^^^ 




Sb(s) pec,LIRG 


56 


48.7 


83.2 


127 


Arp220^ 




S,SyULIRG 


76 


57.0 


96.5 


133 


NGC7252''^ 


Arp 226,The Atoms of Peace 


(R)SA(r) 


63 


25.1 


127.0 


214 


Arp222^ 




SAB(s)apec 


23 


43.4 


66.5 


N/A 


NGC 1700^ 




E4 


54 


90.0 


87.0 


N/A 



-1 
-1 
-1 
-1 
-1 
-1 
-1 




1 
1 
1 
1 
1 
1 
1 

5 
5 
7 
9 
10 
11 
12 
13 



Notes. *^' - a member of the Toomre sequence ( I1977I I; '^' 
compiled sample of angularly-large interacting galaxies 



a member of the X-ray sample of Brassington et al. ( I2007t ; *■'' - a member of our 



Table 2. Parameters of the reduced radio observations at 4.86 GHz. 



Name 



Project ID Amp. cal. Phas. cal. Integr. time 

[h] 



Stp Sp, r.m.s. (TP) r.m.s. (PI) 

[miy] [miy] [juIy/beam] [/ujy/beam] 



164 7+'~'~ 
45 0+"-^ 

1 10-7+0.9 

44 9+0.9 
50.3+'^-^ 



9 q+nr 
^•'^-o.i 

2 r'- 

10 5+"^' 
3.8!«? 

2.7+; 

1.0;^' 

4+o:6 



NGC 660 


AL296 


3C138 


0202H-149 


4.0 


NGC 876/7 


AH314 


3C48 


0202-^149 


5.2 


NGC 2207/IC 2 163 


AH370 


3C286 


0607-157 


5.0 


NGC 5426/7 


AL425 


3C286 


1406-076 


3.9 


NGC 6907/8 


AH314 


3C286 


2128-208 


5.2 


TAFFY 


AC264 


3C48 


0007+171 


1.4 


TAFFY2 


AC592 


3C48/3C138 


0114+483 


5.4 



38.4 
19.2 



-M 



26.4 
16.5 
40.0 
23.3 
32.0 
17.0 
23.0 



13.4 
14.0 
21.2 
14.9 
23.3 
16.0 
19.0 



Table 3. Parameters of the original images from the NRAO VLA Archive* 



Name 



Project ID Frequency VLA config. Integr. time 



[GHz] 



NGC 520 


AY 102 


1.43 


NGC 1700 


NVSS 




NGC 3256 


AS412 


4.86 


NGC 6621/2 


AU23 


1.49 


NGC 7252 


AN090 


4.86 


Arp 220 


AY 102 


1.43 


Arp 222 




1.49* 



c 

D 
BC 
BC 
BC 
C 
D 



Notes. **' from http://www.aoc.nrao.edu/~vlbacald/[ ***' from Condon ( 11987b 



[h] 



1.0 

0.2 
0.1 
1.3 
0.6 



Stp r.m.s. (TP) 

[miy] Luly/beam] 



172.85y 

<5.0 
250.0+i'/ 

28.0!^; 

Q 1+0.9 
'•^-0.2 

334.3!0-5 

3+0.7- 



2.9 



-0 1 



150.0 
350.0 
265.0 
225.0 
34.7 
175.0 
110.0 



of colliding galaxies, they do represent the different stages of 
tidal interaction and merging process, which is suitable for our 
study. 

We reduced the original VLA radio data for the seven inter- 
acting systems listed in Table|2] Some of them were unpublished 
to date, some were published but without any analysis of polar- 
ization data, and for some objects we obtained the final images of 



better quality than those known from the literature (see the next 
section). Using the standard AIPS package and procedures, the 
instrumental polarization was corrected with phase calibrators, 
which were also applied in gain and phase calibration (Table|2]i. 
The flux density scale was established by the amplitude calibra- 
tion sources and the data were self-calibrated in phase, except 
for NGC 5426/7. NGC 6907/8 data were also self-caUbrated in 
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NGC2207;lC2163 VLA 4.86GHz TP+PI B-vect. on DSS blue 



NGC2207/IC2163 VLA 4.86GHz Pl+B %-vect. on DSS blue 




Okpc 



06 16 35 30 25 20 15 

RIGHT ASCENSION (J2000) 

Fig. 1. The total-power contours and B -vectors of polarized in- 
tensity of NGC2207/IC2163 at 4.86 GHz (natural weighting) 
superimposed on the DSS blue image. The contour levels are (- 
3, 3, 5, 8, 12, 20, 35, 80, 150) x 40.0 //Jy/beam. A vector of 10" 
length corresponds to the polarized intensity of 166.7 yuJy/beam. 
The map resolution is 27" x 27" HPBW. 




06 16 35 30 25 20 

RIGHT ASCENSION (J2000) 

Fig. 2. The contours of polarized intensity and B-vectors of 
polarization degree of NGC2207/IC2163 at 4.86 GHz (natural 
weighting) superimposed on the DSS blue image. The contour 
levels are (3, 5, 8, 12, 20, 35) x 21.2 yuJy/beam. A vector of 10" 
length corresponds to the polarization degree of 10%. The map 
resolution is 27" x 27" HPBW. 



amplitude. The combined maps in I, Q, and U Stokes parame- 
ters at both frequencies were then used in constructing the maps 
of polarized intensity, polarization position angle, and position 
angle of observed magnetic vectors (observed E-vectois rotated 
by 90° - the so-called B-vectors). The details of observations, 
obtained r.m.s. sensitivities of the radio maps at 4.86 GHz, and 
total and polarized radio fluxes are presented in Table |2l 

Also at our disposal were reduced high-quality VLA radio 
polarization data at 4.86 GHz for the Antennae galaxies and for 
the Virgo Cluster member NGC 4254 from our previous publi- 
cations (seefor details Chyzy & Beck [2004l Chyzy et al. l2007bl 
Chyzy I2008I I). For the seven other mergers we used the orig- 
inal images from the VLA archive. The details of those data 
and r.m.s. sensitivities in total radio emission are presented in 
Table [3] The short integration time prevented us from calibrat- 
ing the polarized emission for these seven objects. 



3. Results 

In this section we describe radio total and polarized emission 
of our interacting galaxies and compare them with morpholog- 
ical properties in other wavelengths. A detailed presentation of 
resolved systems is followed by concise descriptions of objects 
without well-resolved radio structures. Next, we estimate mag- 
netic field strengths for all the objects investigated. 

3.1. NGC2207/IC2163 

In this pair of spiral galaxies, NGC 2207 appears in the fore- 
ground, obscuring the other one, IC2163. According to the 
numerical simulations, 4 x 10^ yr ago IC2163 passed behind 
NGC 2207 from west to east at a distance of about 20kpc 
(Elmegreen et al. 1995b). Seventeen young stellar clusters were 
discovered in both galaxies (Elmegreen et al. 12001), but the in- 
teraction has not yet led to a prominent starburst. 

The radio emission at 4.86 GHz (Fig. [TJ covers the opti- 
cal extent of the system but displays a hole in the central part 
of NGC 2207 Uke in the map at 1.4 GHz (Condon fT983T l. The 
emission traces the optical western tidal tail of IC2163 and 



peaks at the eastern edge of NGC 2207, where a gravitational 
instability region is suspected, producing a maximum in the in- 
frared emission (Elmegreen et al. I2006I I. In this region, the to- 
tal and polarized emission also have maxima, but not that de- 
gree of polarization, so there is no enhanced production of or- 
dered magnetic field here, and the found radio properties of this 
region could result from a supply of relativistic electrons com- 
ing from the western tidal tail of the IC 2163. This idea is sup- 
ported by accretion of gas from the IC 2163 onto the NGC 2207, 
as suggested by Elmegreen et al. (il995al ). The Hi velocity dis- 
persion is quite high in a large area of the disk of NGC 2207 
(40 - 50kms '). However, we do not see any correspondence 
of regions with enhanced velocity dispersion (cf. Elmegreen et 
al. 1995a) with observed distribution of total or polarized ra- 
dio emission. Therefore, the tidally generated turbulence and/or 
stretching motions along the line of sight required to explain the 
observed pattern of velocity dispersion cannot significantly in- 
fluence the magnetic field strength. 

The regular magnetic field around the central part of 
NGC 2207 have an almost exclusively azimuthal direction (Fig. 
|2]i. Also HI gas forms a similar ring-like structure (Elmegreen 
et al. I1995al l. However, in the southern part of the disk, the B- 
vectors have a larger pitch angle ( 45°) and go out of the disk. 
This may indicate an outflow of magnetized plasma along a tidal 
tail and conversion of random magnetic field to the regular com- 
ponent as in the eastern tail of the Antennae galaxies (Chyzy 
& Beck 2004). With rising distance from the disk, the degree 
of polarization rises to 40%, and the magnetic field becomes 
more regular In this region, the optical and infrared emission 
form long and narrow structures called "streamers" and have the 
same pitch angle as B-vectors (Fig.|2]l. In this part of NGC 2207, 
a large bubble of H i gas extends up to 30kpc out of the disk to 
the south (Elmegreen et al. ll995a] l. 

A similar behaviour of the magnetic field is also observed in 
the eastern part of IC 2163. The degree of polarization becomes 
larger with distance from the disk centre (reaching up to 30%), 
which can indicate stretching of field lines along the tidal tail 
(Fig.EJ. 
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NGC5426/7 VLA 4.86GHz Pl+B %-vect. on DSS blue 
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Fig. 3. The total power contours and B-vectors of polarized in- 
tensity of NGC 5426/NGC 5427 at 4.86 GHz (natural weighting) 
superimposed on the DSS blue image. The contour levels are (- 
3, 3, 5, 8, 12, 20, 35, 80, 150, 250) x 23.3 ^^Jy/beam. A vec- 
tor of 10" length corresponds to the polarized intensity of 83.3 
//Jy/beam. The map resolution is 22" x 22" HPBW. 
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Fig. 4. The contours of polarized intensity and B-vectors of po- 
larization degree of NGC 5426/NGC 5427 at 4.86 GHz (natural 
weighting) superimposed on the DSS blue image. The contour 
levels are (3, 5, 8, 12) x 14.9 //Jy/beam. A vector of 10" length 
corresponds to the polarization degree of 17%. The map resolu- 
tion is 22" X 22" HPBW. 



3.2. NGC 5426/7 

Both galaxies in this binary system are spirals of the same 
Hubble class Sc, with similar sizes (x; 20kpc), and masses 
( 3 X 10'°MO; Fuentes-Carrera et al. .2004) . The evidence for 
weak gravitational interaction are the thin straight strands of the 
outer arms of NGC 5426 that overlap the southern spiral arms of 
NGC 5427 and the straight arm segment in the SW part of the 
spiral structure of NGC 5427 (Fig. [3]l. The velocity field of ion- 
ized gas in NGC 5427 is also slightly distorted (Fuentes-Carrera 
et al. l2004l l. Blackman ( fT982] l found that 38% of the total opti- 
cal flux of the system comes from two adjacent halves of both 
the disks. This suggests obscuration by dust in NGC 5426 that 
situates this galaxy in front of NGC 5427 (Fuentes-Carrera et 
al. 120041 . All these features can be accounted for by an early 
phase of interaction that suggests a closer passage of galaxies in 
5x10' yrs (Fuentes-Carrera et al. l2004l l. 

This scenario of weak tidal interaction fully agrees with the 
radio polarimetric data at 4.86 GHz. The total radio intensity 
closely follows the optical morphology (Fig. [3]). Both galaxies 
have coherent magnetic spiral patterns that closely follow the op- 
tical spiral arms. The absence of polarized emission in between 
the galaxies (Fig. |4| is probably due to beam depolarization of 
two structures of perpendicular B-vectors that overlap there. 

In the most disturbed S W part of NGC 5427 there is no po- 
larized emission, which makes this galaxy slightly asymmetric 



in the polarized intensity (Fig. |4]i. In the SE part of NGC 5427 
the polarized emission occurs outside the density wave. To date, 
such outer magnetic arms have only been observed in the weakly 
interacting Virgo cluster spiral NGC 4254 (Chyzy et al. I2007bl l. 

3.3. NGC 876/7 

The larger component of this pair of galaxies - NGC 877 - is 
a late-type spiral with prominent and slightly disturbed optical 
spiral arms (Fig. |5]). Bright Ha emitting regions are distributed 
along the arms but are not condensed in the central part of the 
disk (Hattori et al. I2004I I. The companion is a small, highly in- 
clined galaxy (78°) at a projected separation of about 50kpc. 

The total radio emission at 4.86 GHz of both the galaxies 
corresponds well to optical and Ha disks and even reveals a 
bridge of low-surface brightness emission (at 180yuJy level) be- 
tween the components (Fig. |5]l. The strong and elongated radio 
structure in the main body of NGC 877 is a galactic bar, from 
which regular magnetic fields stretch out to the N and S. The 
magnetic pattern is very coherent one and resembles those in 
grand-design spiral galaxies. The presence of tidally excited spi- 
ral arms in a system with a small companion is predicted by 
numerical simulations (Toomre'1981). This would explain the 
excess in infrared emission of NGC 877, which led to classify- 
ing this object as LIRG (Sanders et al. I2003I I. The smaller com- 
panion (NGC 876) has more concentrated radio emission, which 
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Fig. 5. The total power contours and B-vectors of polarized in- 
tensity of NGC 876/NGC 877 at 4.86 GHz (natural weighting) 
superimposed on the DSS blue image. The contour levels are 
(-3, 3, 5, 8, 12, 20, 35, 80, 150, 250) x 16.5 yuJy/beam. A vec- 
tor of 10" length corresponds to the polarized intensity of 83.3 
/iJy/beam. The map resolution is 18" x 18" HPBW. 
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Fig. 6. The contours of polarized intensity and B-vectors of po- 
larization degree of NGC 876/NGC 877 at 4.86 GHz (natural 
weighting) superimposed on the DSS blue image. The contour 
levels are (3, 5, 8) x 14.0 /iJy/beam. A vector of 10" length cor- 
responds to the polarization degree of 14%. The map resolution 
is 18" X 18" HPBW. 



may result from high inclination (78°) of this spiral. Surprisingly 
enough, it shows clear hints of the X-shape structure of B- 
vectors, often observed in the edge-on field galaxies (Fig.|5]l. 

The distribution of the polarized emission of NGC 877 is 
asymmetrical (Fig. |6]l. The major part of this emission comes 
from the NW part of the disk, where the spiral arm is also better 
visible in the optical images. The optical morphology in the SE 
part of the disk is only slightly distorted but corresponds to a dis- 
tinct magnetic spiral arm. A blob of polarized emission in the SE 
part of the disk coincides with an appendage in the near-infrared 
emission, which also accounts for asymmetrical appearance of 
this object at this band (Moriondo et al. |1999l l. The companion 
does not appear clearly asymmetrical in the polarized emission, 
indicating a weaker effect from interaction onto this object. The 
weak distortions of magnetic fields also indicate an early type of 
interaction of this system. 

3.4. NGC 6907/8 

NGC 6908 is a small lenticular galaxy superimposed on the NE 
arm of the dominant galaxy NGC 6907. It is most clearly seen 
in near-IR images (Madore et al. 12007!) . The spiral structure of 
NGC 6907 is distorted: its eastern arm is much longer and wraps 
around the entire galaxy (Fig. |7]i. The bulge of NGC 6907 has 
an elliptical shape, and is shifted off the base of the arms. The 
analysis of distortions in the HI velocity field suggests that the 
first collision between galaxies occurred 3.4 + 0.6 x 10^ yr ago, 
leaving behind a triangular HI gas stream to the NE from the 
current position of NGC 6908 (Scarano et al. .2008) . Sanders et 
al. (I2003l l classifies NGC 6907 as LIRG. 

In our maps at 4.86 GHz, the total radio emission is detected 
only from the larger galaxy and concentrates in its central part 
(Fig.EJ. Two protrusions to N and S are likely some background 
objects. The polarized emission (Fig. [8) is more asymmetric and 



dominates outside the galactic bar in the SE part of the disk. 
The polarized data reveal a coherent spiral structure of magnetic 
B-vectors of small pitch angle (0°-30°), even in its eastern part, 
where the optical density waves are weak. In the eastern part of 
the disk the well-ordered magnetic fields seem to be pulled out 
and directed along the tidally stretched spiral arm. 



3.5. The Taffy 

The unusual morphology of this pair of galaxies (UGC 12914 
and UGC 12915), known as the Taffy system, comes from a re- 
cent (x 10^ yr), nearly head-on encounter (Condon et al |1993l l. 
Almost half of their radio emission at 1 .49 GHz have been de- 
tected just from the bridge between the galaxies, which may be 
due to fresh relativistic CR electrons radiating in magnetic fields 
pulled out of the disks. In this system the Hi gas is mostly dis- 
tributed along the bridge (Condon et al. I1993I I, where a large 
amount of molecular gas and dust were also detected (Gao et al. 
r2003;Zhuetal. l2007] l. 

Our re-reduced radio polarization data at 4.86 GHz (Fig.|9l) 
show more details of the magnetic field structure than are pre- 
sented in Condon et al. (1993). We have confirmed that the 
general orientation of magnetic field is along the line join- 
ing the galaxies, which is also a possible direction for passing 
the intruder galaxy UGC 12914. In addition, we found that in 
UGC 12915 the polarized emission is clearly visible on the in- 
teracting side of the disk, whereas the other side is almost com- 
pletely devoid of polarized signal. Either the magnetic field was 
enhanced asymmetrically by the interaction or the magnetized 
ISM on the NE side of the disk was stripped away by ram pres- 
sure during the passage through UGC 12914. The UGC 12914 
shows some patches of polarized signal on the disk side oppo- 
site the bridge, which are associated with ordered magnetic fields 
with large pitch angles of up to 90° (Fig.|9]i. 
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Fig. 7. The total power contours and B-vectors of polarized in- 
tensity of NGC 6907/NGC 6908 at 4.86 GHz (natural weighting) 
superimposed on the DSS blue image. The contour levels are (-3, 
3, 5, 8, 12, 20, 35, 80, 150, 250, 400) x 32.0 yuJy/beam. A vec- 
tor of 10" length corresponds to the polarized intensity of 125.0 
//Jy/beam. The map resolution is 30" x 30" HPBW. 
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Fig. 8. The contours of polarized intensity and B-vectors of po- 
larization degree of NGC 6907/NGC 6908 at 4.86 GHz (natural 
weighting) superimposed on the DSS blue image. The contour 
levels are (3, 5, 8, 12) x 23.3 yuJy/beam. A vector of 10" length 
corresponds to the polarization degree of 8%. The map resolu- 
tion is 30" X 30" HPBW. 



TAFFY(UGC12914/5) VLA 4.86 GHZ TPtPI B-vect. on DSS blue 
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Fig. 9. The total power contours and B-vectors of polarized in- 
tensity of the Taffy (UGC 12914/UGC 12915) at 4.86 GHz su- 
perimposed on the DSS blue image. The contour levels are (-3, 
3, 5, 8, 12, 20, 35, 80, 150, 250, 400) x 17.0pJy/beam. A vec- 
tor of 10" length corresponds to the polarized intensity of 71.4 
/iJy/beam. The map resolution is 16" x 16" HPBW. 



TAFFY2(UGC813/6) VLA 4.86 GHZ TP+PI B-vect. on DSS blue 
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Fig. 10. The total power contours and B-vectors of polarized in- 
tensity of the Taffy2 (UGC 813AJGC 816) at 4.86 GHz superim- 
posed on the DSS blue image. The contour levels are (-3, 3, 5, 8, 
12, 20, 35, 80, 150) x 23.0yL(Jy/beam. A vector of 10" length cor- 
responds to the polarized intensity of 69.4 //Jy/beam. The map 
resolution is 15" x 15" HPBW. 



3.6. The Taffy2 

This pair of post-collision spiral galaxies (UGC 813 and 
UGC 8 16) is similar to the Taffy system and shows a radio bridge 
at 1.40 GHz (Condon et al. 2002). The suspected head-on colli- 
sion was also a recent one that occurred about 5 x 10^ yr ago. We 
call this system the TafFy2. 

Our re-reduced data at 4.86 GHz show that, besides the ra- 
dio bridge, the regular field is discerned in the barred disk of 
UGC 816 and resembles a typical spiral structure (Fig.fTOli. The 
configuration of B-vectors is not as clear in the second galaxy. 



which could be caused by its high inclination (72°). The orienta- 
tion of magnetic field is almost perpendicular to the line joining 
the disks, which is somewhat puzzling, considering the oppo- 
site orientation of magnetic field in the Taffy system. The uncer- 
tainty of magnetic field orientations estimated from the signal- 
to-noise ratio of the polarized intensity is about 10°. On the other 
hand, the foreground rotation measure of the Milky Way is about 
lOradm"^ for the Taffy coordinates and about -50radm"^ for 
the Taffy2 (Johnston-Hollitt et al. 2004). Thus, the uncertainties 
of the magnetic field orientation due to the foreground RM are 2° 
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Fig. 11. a) The total power contours of NGC660 at 4.86 GHz (natural weighting) superimposed on the DSS blue image. A large 
warped ring in the NS direction is visible. The contour levels are (-3, 3, 5, 8, 12, 20, 35, 80, 150, 250, 400, 700, 1500, 4000) x 
26.4 juJy/beam. b) The total power contours and B-vectors of polarized intensity of the disk of NGC 660 superimposed on DSS 
blue image. The contour levels are the same as for Fig. a). A vector of 10" length corresponds to the polarized intensity of 100.0 
yuJy/beam. c) The contours of polarized intensity and B-vectors of polarization degree of the disk of NGC 660 at 4.86 GHz (natural 
weighting) superimposed on the DSS blue image. The contour levels are (3, 5, 8, 12, 20, 35) x 13.4 yuJy/beam. A vector of 10" 
length corresponds to the polarization degree of 1 1%. The maps resolution is 18" x 18" HPBW. 



and 11°. Therefore, neither uncertainties in determination of ori- 
entations of B-vectors nor the foreground RM could cause such 
large differences in the observed patterns of magnetic fields. 



Following from this, magnetic field orientations perpendicu- 
lar to the interaction axis in the TafFy2 shows that this system is 
less taffy-like than the Taffy one. This may indicate less violent 
interaction in Taffy2, hence weaker distortions of magnetic fields 
out of its galactic disks. The rough estimation of kinetic energy 
involved in gravitational interactions are 4.9- 14. Ox 10^^ erg and 
3.2 - 4.6 X 10^^ erg for TaflFy and Taffy2 galaxies, respectively 
(Lisenfeld & VoLk 1201 Oi l. The Taffy2 thus indeed gives a hint 
of weaker interaction. Another possibility for different magnetic 
field orientations is a former collision in Taffy2 giving a longer 
time for galactic rotation to re-establish the global galactic mag- 
netic field (especially in UGC 816, cf. Fig.fTOli. The differences 
in magnetic field configuration may also arise from projection 
effects, provided the magnetic field lines are not straight but are, 
e.g., helical in the bridges (M. Soida, private communication). 
Detailed MHD simulations are needed to evaluate these possi- 
bilities. 



3.7. NGC 660 

NGC 660 is usually recognized as a polar-ring galaxy with a 
large warped ring of gas and stars (Fig. \TTk). It was suggested 
that its morphology results from a central collision and subse- 
quent merging of two spiral galaxies of similar masses that hap- 
pened several Gyrs ago (van Driel et al. ll995l l. However, numer- 
ical simulations can indicate another scenario where the tidal in- 
teractions between two objects without a final merger caused ac- 
cretion of polar-ring material to NGC 660 from the donor galaxy 
(Bournaud & Combes 2003). 

Our map of the total radio intensity at 4.86 GHz (Fig. fTTb) 
displays a smooth emission from a disk with a strong central 
component. The central source has a flux of 85 ^g mJy, which 
agrees with its total flux at 1.41 GHz and the spectral index a = 
0.6 (S oc V"") estimated by Condon et al. (1982). Surprisingly 
enough, the directions of B-vectors form an undistorted 'X-like' 
structure in the outer part of the disk and parallel to the disk in its 
inner part (Fig.fTTt). Such a configuration is typically observed 
in normal late-type spiral galaxies seen edge-on (e.g. Braun et 
al. 2010). We suggest that in this case the second evolutionary 
scenario with tidal accretion, rather than one involving a strong 
disk deformation (expected for the merger), is more likely for 
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this system. The convolution of our maps to the larger beam of 
1' does not reveal any radio protrusion from the polar ring. This 
could be expected, because in either evolutionary scenario, the 
life-time of relativistic electrons at 4.86 GHz (about 10*^ yr) is 
much shorter than the estimated age of the ring (about 10'' yr; 
van Drieletal. [19951 1. 

3.8. Other interacting and merging systems 

NGC4254 — This is a weakly interacting spiral galaxy in the 
the Virgo cluster, which has recently experienced a gravitational 
encounter to have perturbed its spiral arms and triggered a burst 
of star formation (Chyzy et al. l2007bl l. The galaxy shows a po- 
larized emission throughout the disk, particularly strong in its 
southern ridge, shifted downstream of a density wave with up 
to 13yuG as a regular magnetic field component and 20 //G in 
the total field (Chyzy 2008). In the ridge, the dynamo-induced 
magnetic fields are modified by stretching and shearing forces, 
which are likely triggered off by tidal forces, which produced an 
anisotropic component of regular field and enhanced the polar- 
ized emission. 

NGC 4038/9 — This well-known interacting pair of almost 
identical spiral galaxies ("the Antennae") is probably just prior 
to the second encounter, which will lead to their merging in 
2 X 10** yr time (Mihos et al. 1993; Kotarba et al. 2010). The 
polarization properties of this system were studied in detail by 
Chyzy & Beck (12004! see the Introduction). The strongest and 
almost random magnetic fields (30;uG) are located in the overlap 
region of both the galactic disks. In the NE tidal tail, away from 
star-forming regions, magnetic field is highly coherent with a 
strong regular component of 1 jjG, tracing gas shearing motions 
along the tail. The radio spectrum is much steeper in this region, 
indicating ageing of the CR electron population as they propa- 
gate from their sources in the star-forming regions. This system 
was chosen by Toomre dl 9771 1 as the first step in his merging 
sequence. 

NGC 6621 — This galaxy, together with NGC 6622, is in- 
cluded in the Toomre sequence as an intermediate-stage merger 
At 1 .49 GHz we see in our maps only the total power emission 
associated with the optically brightest part of the system. 

NGC 520 — This system represents an intermediate stage 
of galaxy interactions according to Toomre (119771 1. The radio 
emission at 1 .4 GHz is strong in the central part of the system 
and shows a protrusion to the NW along the optically bright 
ridge. According to Stanford & Barcells (119911 ). the system is 
at present 3 x 10^ yr past the encounter and prior to the final 
coalescence. 

NGC 3256 — A kinematic study of this ULIRG system 
(English et al. 120031 ) indicates that it is currently experiencing 
a starburst just prior to the merging. The two tidal tails visible 
in optical images probably formed about 5 x 10^ yr ago and the 
coalescence is estimated to take place in w 2 x 10** yr The ra- 
dio data at 4.86 GHz reveal only total emission peaked at the 
galactic centre, the site of intense infrared and X-ray emission 
(Brassington et al. 120071 ) . 

Arp 220 — This is the closest ULIRG, the merger at the point 
of coalescence (Brassington et al. 120071 ) . The radio emission at 
4.86 and 1 .43GHz is highly concentrated and coincident with the 
optical one. 

NGC 7252 — This well-known merger, called "The Atoms 
for Peace", shows remarkable optical loops and tidal tails sur- 
rounding it (Dopita et al. 2002 ). This is the last system within the 
Toomre (i 19771 ) sequence of interacting galaxies. The star forma- 
tion is strongly concentrated towards the nucleus and has a coun- 



Table 4. Magnetic field strength and regularity for interacting 
galaxies. 



Name 


B/o, 


Br eg 


Field 




L"G] 


^^G] 


Regularity 


NGC 4254 


15 ±4 


7±2 


0.48 + 0.17 


NGC 5426 


11 ±3 


3±1 


0.32 + 0.13 


NGC 5427 


13 ±4 


4 + 2 


0.33 + 0.13 


NGC 876 


11 ±3 


2+ 1 


0.21 + 0.09 


NGC 877 


15 ±5 


4 + 2 


0.24 + 0.09 


NGC 2207 


16 ±5 


6±2 


0.39 + 0.15 


IC2163 


12 ±4 


4±2 


0.33 + 0.15 


NGC 6907 


15 ±4 


3± 1 


0.23 + 0.10 


UGC 12914 


12 ±4 


3± 1 


0.27 + 0.11 


UGC12915 


15 ±4 


2± 1 


0.15 + 0.07 


UGC 813 


13 ±4 


3±1 


0.21 + 0.09 


UGC 816 


15 ±5 


3±1 


0.19 + 0.09 


NGC 660 


16 ±5 


3± 1 


0.18 + 0.08 


NGC 4038 


18 ±6 


4±2 


0.24 + 0.11 


NGC 4039 


12 ±4 


3±2 


0.23 + 0.12 


NGC 6621 


13 ±4 


N/A 


N/A 


NGC 520 


13 ±5 


N/A 


N/A 


NGC 3256 


25 ±8 


N/A 


N/A 


Arp 220 


27 ±7 


N/A 


N/A 


NGC 7252 


12 ±4 


N/A 


N/A 


Arp 222 


5±2 


N/A 


N/A 


NGC 1700 


<6±2 


N/A 


N/A 



terpart in the radio emission at 4.86 GHz. According to Mihos et 
al. ( 119931 ), the SFR had dropped in this object by two thirds from 
the time of coalescence, estimated to occur x 1 Gyr ago. 

Arp 222 — A post-merger galaxy is at a slightly more ad- 
vanced stage of evolution than NGC 7252 (Brassington et al. 
2007 ). The protrusions and plums observed in optical and X-ray 
emission indicate that the system has not relaxed to an ellipti- 
cal galaxy yet. Our analysis of available VLA radio data of this 
system at 1 .49 GHz provides the total emission from the central 
part of the system. 

NGC 1700 — A protoelliptical galaxy possesses a counter- 
rotating stellar core and other symptoms that it formed through 
merging of two spiral galaxies of similar mass (Brassington et al. 
,2007). The best age estimate for the coalescence event is about 
3 Gyr (Brown et al. 2000). No radio emission of this galaxy was 
detected in the NVSS survey, so only the upper limit of radio 
emission at 1.4 GHz could be estimated (Table[3]). 

3.9. Magnetic field strengtii and regularity 

The radio emission detected by us for interacting galaxies is 
composed of nonthermal (synchrotron) and thermal (free-free) 
radiation. The nonthermal emission and its polarized component 
can be used to obtain total and regular (ordered) magnetic field 
strengths (Beck & Krause 2005. see Appendix lAl for details). A 
typical equipartition energy condition between cosmic rays and 
magnetic fields can be applied, especially as spatial resolution of 
the the analysed radio observations are as large as hundreds of 
parsecs. For all galaxies the total and polarized radio intensities 
were integrated over the galactic disks as outlined by the opti- 
cal emission visible in optical B-images. We assumed a typical 
unprojected pathlength through the synchrotron emitting region 
of 1 kpc. However, for NGC 1700, which is perfectly an edge-on 
object (Table [TJ without detectable radio emission, we adopted 
the synchrotron pathlength L of 2 kpc (see AppendixlAl. 
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The nonthermal spectral indexes and thermal fractions 
for NGC 2207/IC 2163, NGC 660, NGC 4038/9, and NGC 4254 
were taken from the separation of radio emission components 
by Niklas et al. ( 1997 ), Chyzy & Beck (l2004l l. and Chyzy et al. 
(12007 b). For other objects, a typical nonthermal spectral index of 
0.8 was used (Niklas et al. II 9971 1. and the radio thermal emission 
was estimated from the FIR data (Surace et al. 2004; Sanders et 
al. 12003 jB rassington et al. 2007; Bushouse et al. 1988; Moshir 
et al. II 9901 1 in the following way. First, we derived global SFRs 
from galactic infrared luminosities at 60 and lOOyum, accord- 
ing to the empirical SFR indicator derived by Kennicutt ( 1998 ). 
They range from less than IMoyr"' (NGC 5426, ARP222, 
NGC1700) up to 15OM0yr-' (ARP220), with uncertainties 
less than 25%. Next, we converted the obtained SFRs to the 
predicted radio thermal emission, according to the modelling 
of free-free radiation by Caplan & Deharveng (1 19861 1. see also 
Chyzy et al. (2007b). The estimated thermal fractions from this 
method are about 0.10 at 4.86 GHz with uncertainties of about 
0.02. Synchrotron emission of galaxies was obtained by subtrac- 
tion of their radio thermal components from the total radio emis- 
sion. 

As a consistency test of the applied approaches to calculating 
galactic magnetic field strengths, we compared the two methods 
of estimating of thermal emission: from the FIR SFR indicator 
and from the separation of thermal and non-thermal radio emis- 
sion. For example for a weakly interacting galaxy NGC 4254 for 
which magnetic properties are well known (Chyzy 120081 ). the 
first method gives £,„, = 14 + 4-fiG, Breg - 6 + 2fiG. From the 
second method we have B,,,, = 15 + 4juG, Breg - T ± 2fiG. 
Therefore, both methods give consistent results that are well 
within estimated uncertainties. We also note that the magnetic 
field strength only weakly depends on the thermal fraction (Eq. 
lA.ll in AppendixlAl). 

The derived strengths of the total and regular magnetic 
fields for all our interacting galaxies are given in Table |4] The 
presented uncertainties include uncertainties of estimated non- 
thermal radio brightness, a factor of 2 uncertainties in adopted 
proton-to-electron density ratio and pathlength L, and the uncer- 
tainty of 0.1 of applied nonthermal spectral indexes. The mean 
values of total magnetic field strength B,o, obtained for the inter- 
acting objects range from 5-6 juG for post-mergers Arp 222 and 
NGC 1700 to 25 - 27juG for ongoing mergers NGC 3256 and 
Arp 220 (Table |4|. This spread may result from the diversity in 
properties of individual galaxies, as well as from the particular 
stages of their gravitational interactions (see the next Section). 
The averaged total magnetic field strength for the whole sam- 
ple is 14 + 5yuG. Excluding the two lowest values for the post- 
mergers (Arp 222 and NGC 1700), which much differ from the 
rest of the sample, the mean strength rises to 15 +4yuG. Both the 
means are more than the typical average total field strength of 
9//G estimated for 74 bright spiral galaxies (Niklas 1995, Beck 
& Krause l2i305l ). As the previous studies of magnetic fields have 
involved all morphological types of galaxies, the cited mean is 
likely increased by including of interacting objects. For instance, 
the average strength of total magnetic field for typical large spi- 
rals in our neighbourhood (the Galaxy and M31) is 5 - 7juG, 
which is significantly lower than for the interacting objects in 
our sample. 

We also calculated the strength of regular magnetic field 
Breg and the field regularity Breg/ Bran, which is a useful mea- 
sure of relative production of regular and random magnetic fields 
(Table 3]). In contrast to the total magnetic field, the strength of 
regular component is similar to those in the non-interacting spi- 
rals, which leads to a mean of 2.5 - 3.5 //G. The enhanced total 



field with relatively normal regular component results in a small 
field regularity (0.27 ± 0.09) in the interacting objects. 

The local magnetic fields are strongest in the centres of 
galaxies and reach 48yuG in Arp 220 and 53yuG in NGC 3256. 
Similar values were also observed in the centres of most star- 
bursting galaxies, as in M 82 (Klein et al. |1988l ). The regular field 
peaks in the regions of weak total radio intensities, far off the 
regions of intense star formation. For example, in the bridges of 
Taffy and Taffy2 galaxies, Breg reaches 1 1 ±3 yuG. This value was 
obtained for the nonthermal spectral index of 1.3 + 0.1 approxi- 
mated from the spectral index maps of Condon et al. ( 119931 ) and 
Condon et al. (2002). Furthermore, we adopted a synchrotron 
pathlength of about 15 kpc corresponding to the width of the ob- 
served radio bridges (cf. Figs. l9l and [TOb. A similar strength of 
regular field was observed at the base of the eastern tidal tail in 
the Antennae system (Chyzy & Beck 2004). Because in all such 
regions the total field is also strong, i.e. 16 ± 4/vG in the bridges 
of the Taffy and Taffy2 systems, the magnetic field can affect 
the local gas dynamics. Actually, in the bridge of the Taffy sys- 
tem, the magnetic energy of 10.2 x 10"'^ erg cm"-* is much larger 
than the energy of thermal gas, which is 1 .4 x 10"'^ erg cm"^ (as 
estimated from the radio thermal emission at 4.86 GHz, assum- 
ing 10** K gas temperature, and the volume filling factor of 0.5). 
The turbulent velocity of Hi gas of about 30kms"' suggested 
by Condon et al. (1993) results in the turbulent energy of ISM 
of 6.6 X 10^'^ergcm"^, again lower than the magnetic energy. 
Although the kinetic energy of H i blobs orbiting the system is 
about two orders of magnitude higher, locally magnetic fields 
can still affect the gas phase, directing the gas motion along the 
magnetic field lines. Similar conclusions are valid for the bridge 
between the Taffy2 galaxies. 

4. Statistical analysis 

It is believed that gravitational interactions between galaxies can 
induce star formation (see Struck 120061 for a review). On the 
other hand, star-forming activity is one of the main factors in de- 
termining dynamo processes and regulating random and regular 
magnetic fields in galaxies (e.g. Chyzy 2008). We have shown 
that the interacting galaxies reveal various radio structures and 
show a range of magnetic field strength (Table |4l). Therefore, 
we investigate how, statistically, magnetic fields and radio mor- 
phologies are connected to the SFR in the interacting objects. 
Having optical sizes of galaxies (corrected for galactic extinc- 
tion and inclination, from HyperLeda), we are able to build sim- 
ilar relations for the surface densities of the SFRs (ESFR). 

In order to relate the properties of radio emission and mag- 
netic fields to different stages of gravitational interactions we 
introduce "the interaction stage parameter", IS, which is sim- 
ply a number indicating the particular stage of interaction and 
galaxy morphological distortion. Toomre ("1977) arranged 1 1 on- 
going mergers according to the rising strength of observed and 
modelled tidal interactions (see Introduction). We assigned them 
consecutive IS values from 1 to 11 (Table [1]). The values 1-9 
correspond to pre-merger stages. IS = 1 describes systems after 
the first galaxy encounter and with significant disk distortions 
(e.g the Antennae and Taffy systems). Galaxies with stronger in- 
teractions, evolving gradually to coalescence, are NGC 662 1/2 
(IS = 5), NGC 520 (IS = 7), and NGC 3256 (IS = 9), all 
displaying strong and concentrated radio emission associated 
with the optically brightest region in the system. The nuclear 
coalescence corresponds to IS = 10 (Arp 220), and the young 
merger remnant to 75 = II (e.g. NGC 7252). We supplemented 
the Toomre sequence with weakly interacting galaxies from our 
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Fig. 12. The evolution of magnetic field strength in interacting 
galaxies. Mean magnetic field strengths are given for the whole 
galaxies (circles), centres (squares), and off-centre regions (as- 
terisks). Arrows indicate upper limits of the field strength. Two 
special stages of interaction are denoted: the first galaxy en- 
counter and the coalescence of the merger nuclei. The error bar 
for NGC 520 is shown. 



sample: we assign IS = -1 to objects before the first encounter 
(e.g. NGC 876/7) and IS = to galaxies close to the (first) 
encounter (e.g. NGC 6907/8) when the tidal tails are launched 
(Sect. |3j Table [U. The other side of the sequence we expand 
with age-advanced mergers but still not relaxed as for Arp 222 
(IS = 12, cf. Brassington et al. l2007] l. Our last object NGC 1700 
is a pre-elliptical system (IS -13) that nevertheless shows ev- 
idence that it was formed from spiral galaxies (Brassington et 
al. (I2007I I. It is clear that the attributed values of IS parame- 
ter may be affected by different interaction geometries or prop- 
erties of individual galaxies. However, to the first approxima- 
tion, they provide a rough indicator of the evolutionary stage 
and "strength" of galaxy interaction. 

4.1. Evolution of magnetic fields 

We explore the evolution of magnetic fields through the interac- 
tion stages of galaxies using the introduced IS parameter (Fig. 
fT2l i. For the weakly interacting galaxies, the strength of mag- 
netic field is almost constant (10- 15 juG). However, we show for 
the first time that there is a distinct trend toward the increasing 
strength of magnetic field as the interaction advances towards 
the nuclear coalescence (up to ISj-iG for Arp 220; IS - 10). A 
completely opposite trend is found at the later stages, i. e. for 
post-mergers, for which the strength of magnetic field decreases 
in a systematic manner. Neither radio emission nor magnetic 
fields are detected for the last galaxy in our sequence, which is a 
proto-elliptical object (NGC 1700). Similar trends for the whole 
galaxies are observed when magnetic field strength is calculated 
only for the central regions of objects (of about 6.5 kpc size) or 
outside of them in the galactic disks (Fig.fTSTi. 

The diverse trend found in magnetic field evolution seems 
to be connected with the star-forming activity in the interacting 
systems. This view is supported by the results of Georgakakis et 



al. (I2000I I. who show for a similar merger sequence that the FIR 
luminosity to molecular hydrogen mass ratio (the star formation 
efficiency) clearly increases close to the nuclear coalescence and 
declines again after merging. The strongly interacting and merg- 
ing galaxies typically have star formation efficiency about an or- 
der of magnitude higher than that of the isolated and weakly in- 
teracting galaxies. Georgakakis et al. also show that, during the 
merging, the enhancement of star formation is moderate within 
the disk, while strong in the nucleus. In accordance with that, 
we observe that the major enhancement of magnetic energy in 
galaxies in our sample also occurs in galactic centres. The field 
strength increases here by more than 170%, whereas within the 
disk the rise is less than 100% (Fig.fTSI). A rise in nuclear lumi- 
nosities has also been found in optical (Laine et al. I2003I I and 
NIR observations with the HST (Rossa et al. l2007l l. 

The H I study of the Toomre merging sequence by Hibbard 
& van Gorkom (119961 ) indicated that the hydrogen gas is pushed 
out from galactic disks into tidal features as the merging process 
advances. Such features are not traced by synchrotron emission 
on our radio continuum maps. Moreover, the large velocity dis- 
persion estimated from H i measurements in weakly interacting 
galaxies (e.g. NGC 2207, Sect. 13. 11 ), likely caused by tidally in- 
duced turbulence or large-scale gas motions, are not associated 
with the observed radio features. Therefore, the star-forming ac- 
tivity must be the main energy source for the radio (synchrotron) 
emission and magnetic fields in merging galaxies. The revealed 
evolution of magnetic fields agrees with the scenario of forma- 
tion of ellipticals from mergers (Toomre & Toomre 1972), be- 
cause there is no detection of synchrotron emission from the el- 
liptical galaxies to date, apart from the AGN phenomena (Sect. 

The evolution of magnetic fields through the interaction 
stages differs significantly from what is observed in the X-rays. 
We do not find any enhancement of magnetic field long after the 
nuclear coalescence, which is observed in the X-ray emission 
in the merging sequence explored by Brassington et al. (2007l l. 
Therefore, it seems that unlike with X-rays, there is no process 
connected to Type I supernova that could re-create the magnetic 
field in the merger remnants, as was proposed to account for the 
X-ray observations (Sect.[T]i. Magnetic fields and cosmic rays re- 
quire powerful Type II supernova explosions for their generation 
(Beck et al. 1 1 9961) . Therefore, our findings suggest that the main 
production of magnetic fields in colliding galaxies must termi- 
nate close to the coalescence, after which magnetic field diffuses 
or is kept at the level of the turbulent energy of ISM. In either 
case, for the radio frequencies considered in this paper, the syn- 
chrotron emission should fade in proto-elliptical (merged) galax- 
ies after about 10*^ yr if there are no sources of CR energization. 
The fading synchrotron emission in post-merger systems is too 
weak to be detected using the currently available radiointerfer- 
ometric data. New radio instrumentation operating at extremely 
low frequencies such as LOFAR could possibly detect the faint 
synchrotron emission from the old population of CR electrons, 
revealing such "relic" magnetic fields. 

4.2. Magnetic fields versus star formation rate 

In Fig. [T3h we show total magnetic field strength related to the 
ESFR for 22 galaxies for which radio emission could be de- 
termined (except for NGC 6908 and NGC 6622, see Sect. [3]). 
We also plotted mean magnetic field strengths in the galac- 
tic centres (Fig. \T3h ) and off-centre regions (Fig. [T3b). We 
found that the statistical relation for the entire galaxies is weak 
(the Spearman correlation coefficient p = 0.49) with an in- 



11 



R. Drzazga et al.: Magnetic field evolution in interacting galaxies 




0.0001 0.001 

LSFR[Msun/kpc"2/yeai-] 

Fig. 13. Total magnetic field strengths (circles) versus local SFRs 
for the whole interacting galaxies a), their centres b), and off- 
centre regions c). The random and regular components of mag- 
netic fields are also provided (squares and asterisks, respec- 
tively). An expected evolutionary track of interacting systems 
is also indicated. 



dex of fitted power law a = 0.20 + 0.02. The dependence is 
stronger in galactic centres (0.27 +0.03) and weaker outside them 
(0.14 + 0.03). We notice that, within a wide range of the ISFR 
of 10"^ - 10"^ Mq yr"' kpc"^, the magnetic field strength is ac- 
tually spread over 10 - 20yuG without any discernible tendency. 
In fact, the determined slopes are established mainly thanks to 
the galaxies of weakest and of strongest SSFR, which are all 
members of particularly strongly interacting systems (IS > 9). 

We also constructed diagrams similar to Fig. [T3] but with 
the global value of SFR. The slopes and correlations are sta- 
tistically similar (e.g. a = 0.21 + 0.02, p = 0.63 for the en- 
tire galaxies). This therefore seems to be a general phenomenon 
that magnetic fields are only weakly dependent on galactic star- 
forming activity, showing a wide spread of estimated strengths. 
The wider the range of considered SFRs, the better defined and 
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Fig. 14. Asymmetry parameters for total {At', squares) and polar- 
ized radio emission (circles; Ap) versus global SFR for interact- 
ing galaxies. The values of offset angles [^t]( 4>p) for asymme- 
try axes in total and polarized emission are given in parentheses. 
Typical uncertainties of asymmetry parameters are less than 5%. 



stronger the influence on magnetic fields. This result can ex- 
plain the former findings of Hummel ( 1981) that the radio con- 
tinuum properties of interacting galaxies differ from those for 
the isolated objects only in the central parts of galaxies, and not 
in their disks. His analysis ignores, however, all the advanced 
mergers that we included in our sample (NGC 3256, Arp 220, 
Arp222, NGC 1700). Similarly, as shown by Kennicutt et al. 
(1987), only particularly strongly interacting galaxies show a 
distinct enhancement in star formation, as appraised from the 
Hot and FIR emission. 

No mutual relation is found between B and SFR (Fig.[T3ll. as 
observed for B and IS (Fig.[T2li. It is because both the advanced 
mergers and the weakly interacting galaxies have low SFRs, and 
they are undistinguishable in Fig. [13] Only when the magnetic 
field is set against the stage of interaction (Fig. \T2[ . is the full 
evolution of B revealed. We plot in Fig.[T3]an evolutionary track 
of exemplary interacting system as based on our results. 

There is a hint of negative correlation between the regular 
component of magnetic field and the ESFR for entire galaxies, 
p = -0.55, and the power-law fit is -0.25 + 0.07 (Fig.[nii)- The 
same trend was more distinct for the magnetic field regularity of 
individual regions of the weakly interacting galaxy NGC 4254 
(Chyzy et al. I2007bl l. Magnetic field tangling in active star- 
forming regions and production of random magnetic fields by 
the local-scale dynamo process (Brandenburg & Subramanian 
2005) result in a lower field regularity. This may also be the case 
over the galactic scales in the analysed interacting objects. 

4.3. Asymmetries 

Gravitational interactions are able to severely distort galactic 
disks (Sect. [U. Here, we investigate how they modify distribu- 
tions of radio total and polarized emission at 4.86 GHz. First, 
we divided each galaxy into two parts by defining an axis in the 
N-S direction that crosses the galactic centre. We integrated the 
total radio intensities over the obtained galactic halves (Ti, T2) 
and determined the ratio aj = TxliJi + T2), with a convention 
of Or > 0.5. Next we rotated the radio map by 30° around the 
galactic optical centre and determined 07- once again. The whole 
procedure was repeated with the step in the rotation of 30°, until 
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Fig. 15. Asymmetry parameters for total (At', squares) and po- 
larized radio emission (circles; Ap) versus the interaction stage 
IS . Typical uncertainties of asymmetry parameters are less than 

5%. 



ual galaxies may blur the effect of possible enhancement of the 
SFR or 2SFR by tidal interactions. 

However, asymmetries are different before and after the first 
galactic encounter (Fig. [TsT l. For all the weakly interacting ob- 
jects (IS = -1; Table [1]), the asymmetries are small with the 
mean values Aj = 0.60 ± 0.03 and Ap = 0.65 + 0.07. For all the 
objects after the first encounter (IS = -Hi), the asymmetries are 
systematically stronger with their means At = 0.68 + 0.03 and 
Ap = 0.87+0. 11 . In polarization they attain the highest measured 
values {Ap > 0.85 for Taffy, Tafry2, and NGC4038). The excep- 
tion from this behaviour is the galaxy NGC6907 (IS = 0; Fig. 
[TsT l. which is close to the moment of collision with NGC 6908 
(Sect.[3]l. However, otherwise than in the other galaxy pairs, this 
system is likely a minor merger and apparently the low-mass in- 
truder (NGC 6908) is not able to distort a much more massive 
galaxy. We thus conclude that the asymmetries visible in the ra- 
dio emission in general only weakly depend on the global or 
local SFRs, while they increase (especially in the polarization) 
with advancing of interaction. As a result, the asymmetry in the 
radio polarized emission could be yet another indicator of an on- 
going merging process. 



the map involved all orientations (six positions). Finally, we de- 
termined the galaxy radio asymmetry At as the highest of all qt 
values. We also calculated the difference ^^ between the origi- 
nal position angle of the major axis of the galaxy (Table [U and 
the position angle of the major axis after rotation, which gave 
the asymmetry Ap- A similar procedure was applied to the radio 
polarization map, which likewise yielded the polarization asym- 
metry Ap and the offset angle of the polarization asymmetry axis 
(pp. To estimate the accuracy of obtained values of asymmetry 
parameters, we rotated the axes which gave asymmetries At and 
Ap by ±15°. The uncertainties approximated as deviations from 
the original asymmetry parameters are small, typically less than 
5%. All the above calculations were performed for those inter- 
acting galaxies for which total and polarized radio emissions 
were resolved and not too much affected by overlapping with 
a companion galaxy (Fig.fT4li. 

The obtained wide range of values of offset angles (cpp, (f>p, 
presented in Fig.[T4b indicates that the asymmetries in total and 
polarized emission are not correlated with the apparent orien- 
tations of galactic disks. However, offset angles are within 60° 
from each other The observed asymmetries are therefore not 
caused by projection effects and must be physically connected 
with galaxy morphologies as modified by the interaction pro- 
cess. 

Among the analysed galaxies there is a trend that the asym- 
metry in distribution of polarized emission (mean of 0.75) is 
stronger than in the total radio intensity (mean of 0.64; Fig.[T4l). 
The ordered magnetic field is therefore much more sensitive to 
morphological distortions than the total field. Even in the weakly 
interacting galaxies (e.g. NGC 876/7), the polarized emission, as 
well as the orientation of magnetic field B-vectors (cf. Sect. [3]), 
are sensitive to global galactic distortions and anomalous gas 
flows. We also notice that the higher the distortion in the total 
emission (Ap), the greater is the difference between the asym- 
metry in total and polarized emission (Ap). 

There is no obvious dependence of the strength of asymme- 
tries on global SFRs for galaxies (p = -0.35; -0.15 for total and 
polarized radio emission, see also Fig.[T4l). We also constructed 
a diagram similar to Fig. [T4]but with the ESFR instead of global 
SFR. Also in this case the results are similar (p = -0.03, -1-0.12). 
It seems that differences in the gas content and mass of individ- 



4.4. Radio-infrared relation 

The process of tidal interaction and induced star formation may 
cause an enhancement in the infrared luminosity. This may lead 
to formation of strong infrared emitters as LIRG or ULIRG ob- 
jects that are believed to be associated with colliding galaxies 
(Sect. [TJ. In order to investigate whether radio emission also 
is correspondingly enhanced in such cases, we constructed a 
radio-FIR relation for our interacting systems and compared it to 
the one determined for a sample of normal galaxies of different 
Hubble types (Fig.fTSTl. We used the radio surface brightness of 
interacting objects estimated from our maps at 4.86 GHz and the 
infrared surface brightness at 60 /vm derived from Surace et al. 
(120041), Sand ers et al. (2003;, Bushouse et al. (fT988] l. and Moshir 
et al. d 1 9901 1. For galaxies with only 1.4 GHz flux available we 
scaled it to 4.86 GHz assuming a typical spectral index of 0.8. 
We did not include two of the weakest radio objects Arp 222 
and NGC 1700 for which 60fj.m infrared data are unavailable. 
Typical uncertainties of obtained values are indicated in Fig. [16] 
and were derived from uncertainties of radio fluxes (Tables|2]and 
[3]l and 60yum fluxes (given in the mentioned FIR catalogues). 
For the sample of normal galaxies radio data at 4.85 GHz were 
taken from Gioia et al. ( 1982) and Chyzy et al. ( I2007al ). whereas 
60jum data came from Helou & Walker d 19951 1 and Moshir et al. 
((79901) . 

Our orthogonal fit to the radio-FIR relation for the sample 
of normal galaxies gives a power-law index of 0.98 + 0.04. The 
interacting galaxies follow this relation (Fig.fTSIl and are located 
in the group of the strongest radio and FIR emitters. We note, 
however, that the brightest galaxy in our sample (Arp 220) is 
shifted slightly down from the main trend indicating possible 
flattening of the relation for relatively young mergers. A more 
numerous sample is needed to confirm this hint. Recently we 
have shown that low-mass dwarf galaxies, which are by orders 
of magnitude weaker optical and radio emitters than the inter- 
acting objects, also follow a similar power-law relation, with a 
slope of 0.91 + 0.08 obtained for 2.64 GHz radio data (Chyzy et 
al. l20m i. 

In fact, the tidally induced intense star formation and mor- 
phological distortion in the studied interacting galaxies do not 
seem to change processes responsible for radio-FIR correlation. 
This is quite interesting, as some of those systems are in the 
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process of evolving and transformation towards SO or elliptical 
merger remnants. The large-scale shocks, galactic large-scale, 
or small-scale dynamos, and processes mutually transforming 
magnetic field components (e.g. by compression and shearing, 
Sect. l4.2b therefore result in the production of magnetic energy 
and cosmic rays that are still well-balanced by production of 
thermal energy and dust radiation. 

Interestingly enough, according to our study, the deviations 
of the Taffy and Taffy2 galaxies from the main trend of the radio- 
FIR relation at 4.86 GHz are small, 0.3 dex, and comparable to 
other objects from our sample (Fig.fTSIl. However, according to 
Condon et al. (119931 12002), these galaxies are twice as radio- 
loud as they ought to be if following the radio-FIR relation con- 
structed at 1.4 GHz. The difference between 4.86 and 1.4 GHz 
is probably due to a more extended emission in the bridge at 
1.4 GHz and the steep radio spectral index (up to 1.3-1.4). The 
reason for such strong radio emission at 1.4 GHz is not clear 
According to Condon, the observed radio excess may come from 
the CRs electrons leaked from the galactic disks and trapped 
in enhanced (stretched) magnetic fields in the bridge between 
galaxies. 

However, Lisenfeld & Voelk (120101 1 argue that in the Taffy 
systems, the kinetic energy of colliding gas probably brought 
about large-scale shock waves that accelerated CRs, as well as 
enhanced radio emission in the bridges. As this is not related 
to star formation, a deviation from the radio-FIR relation results. 
They assume in their modelling that the shocks in the Taffy disks 
correspond to the one in the middle-age supernova remnants and 
therefore could be described by the Sedov phase. However, ac- 
cording to our analysis of polarization data in the Taffy system 
(Sect. [U, the magnetic field in the bridge is well-ordered and 
aligned with the direction of gas flows between disks. Thus the 
field configuration clearly resembles a very young supernova 
with magnetic field vectors aligned with gaseous outflows, well 
before the Sedov phase. As a result, the configuration of ordered 
magnetic fields indicates that the physics involved in explaining 
the Taffy bridge might not be applicable here. Alternatively, we 
deal here with some other processes of field ordering and com- 
pression, which work differently than in the supernova remnant. 
Sensitive observations of gamma rays could independently as- 
sess the CR population, indicate their origin, and account for the 
strong radio emission in the bridge. In the Taffy2 system situ- 
ation is different because magnetic field is not aligned with the 
bridge (Fig.fTOb. This may indicate a better correspondence with 
the Lisenfeld & Voelk model. 

4.5. The impact of merger's magnetic fields 

Gravitational forces could play an important role in generation 
of magnetic field in the early-universe. If seed magnetic fields 
were produced by, e.g., Biermann battery or Weibel instability 
then small-scale dynamo during gravitational collapse of galax- 
ies and accretion driven turbulence could strongly enhance the 
seed field (Schleicher et al. .2010) . Further evolution of galaxies 
could build up stronger, galactic-scale magnetic fields by small 
and large scale dynamos as discussed by Arshakian et al. (I2009I I. 
Gravitational forces during galaxy interactions and galaxy merg- 
ers could also play an important role in the magnetization of 
IGM. Although the strongly interacting objects represent just a 
small fraction of galaxies in the local Universe, < 300 interact- 
ing objects within 50 Mpc distance according to the Vorontsov- 
Velyaminov et al. (2001) catalogue, they may have dominated 
the epoch of massive galaxy formation in the early Universe 
(Conselice l2006l l. Such galaxies are able to redistribute magnetic 
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Fig. 16. The radio-FIR correlation diagram for interacting 
(squares) and normal galaxies of different Hubble types (aster- 
isks). The surface brightness at 4.86 GHz and at 60/vm (Jy/D') 
are used (see text for details). The solid line is an orthogonal fit 
to normal galaxies, and dashed lines represent simple X vs. Y 
and Y vs. X regressions. Error bars for interacting galaxies are 
marked. 



fields in IGM along tidal tails or bridges, as we in fact observe, 
e.g., in the Antennae and Taffy galaxies (Sect. [3]l. The lower 
limit on IGM magnetic fields of 5 x 10"'^ yuG was constrained by 
Tavecchio et al. (.2010) based on analysis of TeV and GeV emis- 
sion from blazars. Recent analysis of FERMI-LAT observations 
of blazars and the absence of their TEV envelopes can suggest 
that extragalactic magnetic fields must fill IGM with a large vol- 
ume factor of * 60% (Dolag et al. 120101) . If such magnetic fields 
are not primordial, then they must have been redistributed from 
local sources, like colliding galaxies. 

To address the idea of magnetization of IGM by interact- 
ing objects, we approximate the zone around these galaxies that 
might be filled with magnetic fields simply by the largest extent 
of their observed H i envelopes. The neutral hydrogen of collid- 
ing galaxies often reveals large protrusions that are not optically 
visible (Hibbard & van Gorkom 1996). If such gas is also accom- 
panied by its partly ionized component, then the frozen magnetic 
field can be advected with gas flows and brought out from galac- 
tic disks. According to the collected literature data (Table [TJ, 
we find that such outflows are up to about 200 kpc in size for 
the most disrupted objects in our sample (e.g. Arp220). This 
value is compatible with the results of numerical simulations of 
high-redshift mergers (Bournaud et al. 1201 Oi l. The tidal inter- 
actions (estimated to operate within about 200 kpc) are there- 
fore not sufficient to seed vast volumes of IGM with magnetic 
fields, as required by the Dolag et al. conclusion. Either unknown 
transport processes redistribute magnetic fields of distant merg- 
ers out of expected H i envelopes, or the proposed magnetization 
of space has a different origin. Even so, interactions can effi- 
ciently magnetize the merger surroundings. This could explain 
relatively strong magnetic fields (» lOj-iG) found around host 
galaxies of quasars with Mg II absorption systems up to redshift 
z = 1.3 (Bernet et al. 2008). Chyzy et al. (12011!) suggest that 
massive galaxies in the redshift z=3-5 could have magnetized 
IGM by supernova explosions and galactic winds to the extent of 
~ 200 - 500 kpc. The process of strong gravitational interactions 
and associated magnetized outflows thus exerts a similar impact 
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Table 5. Deflection angle of UHECR in a colliding galaxy sys- 
tem. 
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to properly associate observed UHECRs with the sites of their 
origin. 

We note that the magnetized halo around mergers, however 
large in extent, would only marginally impact the UHECR prop- 
agation (deflection up to about 1°), owing to dispersed, weak 
magnetic fields (Table |5j. The propagation of protons through 
the Galaxy would also result in a small deflection of up to about 

r (for Bran ^ 4 fiG). 

5. Summary and conclusions 

We studied the influence of gravitational interactions on the dis- 
tribution of radio emission and the properties of magnetic fields 
in galaxies. We perused archival radio data from the VLA at 
4.86 GHz and 1 .4 GHz for a sample of 24 galaxies in 16 interact- 
ing systems. Different stages of interaction were approximately 
described by the interaction stage (IS) parameter (Sect. |4]l. The 
main conclusions are as follows. 



on IGM. New observational methods, such as a rotation mea- 
sure synthesis, applied to background sources shining through 
merger remnants, could reveal such magnetized cocoons around 
mergers. New instrumentation such as the EVLA, LOFAR, and 
SKA could be used for such detection experiments to confirm 
these predictions. 

A vital issue in high-energy astrophysics is to identify sites 
of origin of UHECRs. There may be problems arising from de- 
flection of UHECR trajectories by intergalactic magnetic fields 
(Yoshiguchi et al. 2003; Neronov & Semikoz 2009; Ryu et al. 
120 lOj . Most interacting galaxies from our sample are within 
lOOMpc distance and their spread-out magnetic fields could 
in principle afi'ect propagation of UHECR on their path to the 
Earth. To estimate a deflection angle S between the incoming 
direction of UHECR and the sky position of their sources, we 
follow the approach of Neronov & Semikoz ( .2009 ), assuming 
that UHECR protons with energy of 10^" eV propagate through a 
magnetized volume of pathlength Lb- In the case of regular field, 
its strength B^eg influences the deflection angle, while in the case 
of turbulent field, its strength £,„„ as well as its coherence length 
Lbc have to be known. We divide our estimations into those for 
disks of interacting galaxies, for strongly magnetized outflows 
(as bridges or tidal tails), and for envelopes of merging galax- 
ies filled with Hi gas and spread out magnetic fields. We take 
Lbc - 50 pc for galactic disks (a commonly adopted value) and 
the scaled values for the bridges and haloes (Table|5]l. 

We obtained the largest deflections of 23° for the protrusions 
from colliding galaxies filled with the regular field such as in 
the Taffy system (Table |5]). In the case of galactic disks with, 
e.g., regular magnetic fields stretched along spiral arms (like in 
the weakly interacting NGC2207, Sect. [3]l the deflection angle 
is also large (16°). The more disrupted (random) fields are here 
less important (6 ~ 3°). The probability of an individual UHECR 
particle to cross the merger and deflect its way to the Earth is 
very low, < 10"^, taking all 300 interacting galaxies into ac- 
count within 50 Mpc distance, and 200 kpc size of their magne- 
tized haloes. However, large-scale shocks in colliding galaxies 
and centres of merging systems were also proposed as possi- 
ble sources of UHECRs ( Bhatta charjee & Sid llOOOl Giller at 
al. | 20031 Biermann et al. I2009I I. If interacting galaxies gener- 
ate some UHECRs by large-scale shocks or AGNs within them, 
their magnetic fields would severely deflect their trajectories (up 
to a: 23°). Propagation of UHECRs through the disk or magne- 
tized outflows (bridges, tidal tails) would make it very difficult 



- The estimated total magnetic field strengths are from 5 f^G 
(ARP222) to 25-27 yuG (NGC3256, ARP220), with a mean 
of 14 + 5 fiG, greater than for bright non-interacting galaxies. 
The mean field regularity of 0.27 + 0.09 is less than for typ- 
ical spirals, which indicates enhanced production of random 
magnetic fields in the interacting objects. 

- Detailed analysis of NGC2207/IC2163 reveals a region of 
enhanced radio total and polarized emission in the eastern 
part of NGC2207, which can suggest a supply of CR elec- 
trons originating in IC 2163. Our studies of polarized prop- 
erties of this system contradict the earlier hypothesis that the 
enhanced radio emission in this region is due to compression 
of magnetic fields (Elmegreen et al. ll995al) . 

- The structure of magnetic field B-vectors in the polar-ring 
galaxy NGC 660 resembles an X-shape pattern, typical of 
late-type edge-on galaxies, which challenges some earlier 
suggestions that this is a post-merger system. 

- Even in weakly interacting galaxies, the polarized emission 
and the orientation of magnetic field B-vectors are very sen- 
sitive tools for revealing global galactic distortions (Taffy 
galaxies) or the direction of gas flow in tidally stretched spi- 
ral arms and tidal tails (NGC 2207, NCG 6907, NGC 4038). 

- We built a tentative scenario of magnetic field evolution in 
interacting galaxies. For weak interactions the strength of 
magnetic field is almost constant (10 - 15juG), then it in- 
creases up to 2x with advancing of interaction to later stages, 
and decreases again, down to the limiting value of 5 - 6 fiG, 
for old post-mergers. We suggest that the main production of 
magnetic fields in colliding galaxies must terminate close to 
the nuclear coalescence, after which magnetic field diffuses 
or is kept at the level of turbulent energy of ISM. The re- 
vealed evolution of magnetic fields agrees with the scenario 
of formation of ellipticals from mergers (Toomre & Toomre 
1972), 

- Magnetic field strength is weakly regulated by the galactic 
star-forming activity (e.g. SFR), showing a wide spread of 
estimated strength values (with a slope a = 0.21 + 0.02 and 
correlation p - 0.63). The dependence is stronger in galac- 
tic centres and weaker outside. There is also a hint that the 
regular field is anticorrelated with SFR (p = -0.55). 

- Asymmetry in distribution of polarized emission is stronger 
(Ap - 0.75) than in total radio intensity (At - 0.64), which 
can indicate that the regular component of magnetic fields 
is much more sensitive to morphological distortions induced 
by tidal interactions than the random magnetic field. This 
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means that, the asymmetry in the radio polarized emission 
could be yet another indicator of an ongoing merging pro- 
cess. The asymmetries visible in the radio emission do not 
statistically depend on the global or local SFRs, while they 
increase (especially in the polarization) with advancing in- 
teraction. 

- The constructed radio-FIR relation for interacting and non- 
interacting galaxies shows that the tidally induced intense 
star formation, morphological distortions, and processes mu- 
tually transforming magnetic field components (e.g. by com- 
pression and shearing) still lead to producing magnetic en- 
ergy and cosmic rays that are balanced by the production of 
thermal energy and dust radiation. The interacting galaxies 
constitute the group of statistically strongest radio and FIR 
emitters. We also argued that the configuration of ordered 
magnetic fields in the bridge of the Taffy system indicates 
that the modelling of Lisenfeld & Voelk (2010), adopted to 
account for the strong radio emission in the bridge, might not 
be applicable here. 

- The process of strong gravitational interactions can effi- 
ciently magnetize the merger surroundings (likely up to 
about 200 kpc), exerting a similar impact on IGM to its 
magnetization by supernova explosions and galactic winds. 
However, tidal interactions are likely insufficient to fill IGM 
with magnetic fields with a large volume filling factor of 
X 60% as recently suggested by Dolag et al. (2010). If in- 
teracting galaxies generate some UHECRs by large-scale 
shocks or AGNs, their propagation through the disk or mag- 
netized outflows (bridges or tidal tails) can deflect them up 
to 23° and make association of observed UHECRs with the 
sites of their origin less obvious. 

Among the interacting systems, those at the final stages of 
merging as well as the post-merging systems (as elliptical galax- 
ies) are the least studied ones with respect to magnetic field evo- 
lution. After the interaction-induced starburst declines, the pro- 
duction of cosmic ray electrons ceases and cannot further en- 
lighten magnetic fields (Sect 14. It . Such fields may still exist in 
the mergers, before they disperse because of turbulent diffusiv- 
ity, but could only be discovered by the Faraday rotation meth- 
ods applied to the background sources. Thus there is a need for 
very high-resolution and sensitive radio spectropolarimetric data 
from interferometers, not yet available. Future EVLA, LOFAR 
and SKA instrumentation could reveal for the first time such fad- 
ing or relic magnetic fields, if they are actually present in the 
post-merging systems. 
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Appendix A: Derivation of magnetic field strengtfis 

The strength of the total B,ot and regular Breg magnetic field can 
be estimated using radio polarimetric observations from the total 
synchrotron intensity /„ and its degree of linear polarization P„,/, . 
After taking the equipartition between the energy densities of the 
magnetic field and cosmic rays into account, the total magnetic 
field is given by (Beck & Krause 2005 ) 
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where Kq is the ratio of proton to electron number densities, a„ 
the mean synchrotron spectral index, L the pathlength through 
the synchrotron emitting medium, Ep the rest energy of the pro- 
ton, and 
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Constants C2 and cj, are geometrical corrections that depend on 
magnetic field orientation. For a region where the field is to- 
tally regular and has a constant inclination / with respect to the 
sky plane, cj, = (cos/)""^'. For a completely turbulent field, 
C3 = (2/3)""^'*^^. If the synchrotron intensity is averaged over 
a large volume, such as the whole galaxy, the value of c^ has 
to be replaced by its average over all occurring values of /. For 
instance, the strength of the mean regular magnetic field in the 
galactic disk can be obtained from the estimated nonthermal de- 
gree of polarization (Segalovitz et al. |1976b : 
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qiKi+a„) - Breg/B,urb, 6 is the azimuthal angle, and B„ub the 
turbulent component of the magnetic field. 

We assume Kq »; 100 in Eq. lA.li which is consistent with 
acceleration of cosmic rays in supernova remnants and with the 
local Galactic cosmic ray data. For synchrotron sources, which 
have a volume filling factor the pathlength through the source L 
has to be replaced by L x /. However, the equipartition magnetic 
field strength only weakly depends on / e.g. its decrease from 1 
to 0.8 results in the increase of field strength by less than 6%. 
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